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ABSTRACT 

We present 5.2-37.2 fim spectroscopy of the edge-on circumstellar disk CRBR 
2422.8-3423 obtained using the InfraRed Spectrograph (IRS) of the Spitzer Space 
Telescope. The IRS spectrum is combined with ground-based 3-5 /xm spec- 
troscopy to obtain a complete inventory of solid state material present along the 
line of sight toward the source. Archival JHKg imaging as well as 350 iim CSO 
mapping, 850 //m SCUBA mapping and 3 mm OVRO interferometry is used to 
obtain a set of spectro-photometric data covering 1.2-3000 /xm. The ices observed 
toward CRBR 2422.8-3423 are compared to archival ISOCAM-CVF 5-16 /xm ice 
spectra of other nearby sources within 2'. We model the object with a 2D ax- 
isymmetric (effectively 3D) Monte Carlo radiative transfer code using all the 
available observations to constrain the source geometry and dust composition. 
In particular, the location of the observed ices in the disk and envelope material 
is included in the model. It is found that the model disk, assuming a standard 
flaring structure, is too warm to contain the very large observed column density 
of pure CO ice, but is possibly responsible for up to 50% of the water, CO2 and 
minor ice species. In particular the 6.85 /xm band, tentatively due to NHj , ex- 
hibits a prominent red wing, indicating a significant contribution from warm ice 
in the disk. The shape of the CO2 bending mode suggests an interaction with 
up to 20% of the CO ice. It is argued that the pure CO ice is located in the 
dense core Oph-F in front of the source seen in the submillimeter imaging, with 
the CO gas in the core highly depleted. Up to 50% of the CO ice embedded 
in water or CO2 ice (no more than 20% of the total amount of CO) may still 
be located in the disk, assuming constant abundances of these types of CO ice 
throughout the system. Discrepancies among the strength of different water ice 
bands are discussed. Specifically, the observed water ice libration band located 
at 11-13 /im is significantly weaker than that of the model. The model is used to 

^SRON National Institute for Space Research, PO Box 800, 9700 AV Groningen, The Netherlands 
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predict which circumstances are most favourable for direct observations of ices in 
edge-on circumstellar disks. Ice bands will in general be deepest for inclinations 
similar to the disk opening angle, i.e. ~ 70°, except for very tenuous disks. Due 
to the high optical depths of typical disk mid-planes, ice absorption bands will 
often probe warmer ice located in the upper layers of nearly edge-on disks. The 
ratios between different ice bands are found to vary by up to an order of mag- 
nitude depending on disk inclination due to radiative transfer effects caused by 
the 2D structure of the disk. Ratios between ice bands of the same species can 
therefore be used to constrain the location of the ices in a circumstellar disk. 

Subject headings: Radiative transfer — circumstellar matter — stars: individual 
(CRBR 2422.8-3423) — ISM: molecules — infrared: ISM 

1. Introduction 

It has been realized in the last decade that circumstellar disks surrounding young low- 
mass stars are a common occurrence. Along with the discovery of a large number of extra- 
solar planets orbiting older stars, significant observational evidence indicates that planets 
and planetesimals are forming in circumstellar disks within 10^ years. Understanding the 
chemical structure and evolution of the molecular material of circumstellar disks is crucial 
because it determines the initial conditions for the composition of planets. This is especially 
important for the study of icy bodies in the Solar system, such as comets and Kuiper belt 
objects. 

Many recent modeling and observational results indicate that the abundances of molec- 
ular gases in circumstellar disks around young low-mass stars arc significantly lower than 
in general molecular clouds, including species such as CO, HCN and HCO+ (Dutrey et al. 
1997; van Zadelhoff et al. 2001; Thi et al. 2004). Although both photodissociation in the 
surface layer and freeze-out onto dust grains may contribute to the depletion, it is believed 
that freeze-out dominates in the disk mid-plane. Some chemical models also suggest that a 
unique ice chemistry, different from that of dark molecular clouds, may exist in circumstellar 
disks (Aikawa et al. 1997). In particular, these authors predict a large abundance of CO2 
ice in typical circumstellar disks. 

Directly observing ices in circumstellar disks is challenging. The primary reason is that 
cold ices can only be observed via absorption bands in the mid-infrared wavelength regime 
and therefore require the presence of a bright infrared source behind the icy material. The 
best candidates for observations of ices in disks are circumstellar disks viewed close to edge- 
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on, where the sight-hne toward the inner bright regions passes through the disk mid-plane, 
which presumably is abundant in ice. A possible exception to this may be a binary young 
star, in which the sight-line toward one component passes through a disk surrounding the 
other. This may be the case for the T Tau binary (Beck et al. 2001; Hogerheijde et al. 
1997), but binarity in this scenario may significantly complicate the structure of disks in the 
system through dynamical and radiative interaction. Circumstellar disks viewed at very high 
inclination (> 75°) become very faint at mid-infrared wavelengths due to the large optical 
depth through the disk mid-plane. The effect this may have on observable ice absorption 
bands is not clear, but observations of ices in this case may not probe deeply into the 
disk mid-plane. The sensitivities of ground-based 8m class telescopes and the Spitzer Space 
Telescope are now high enough to directly observe ices in a range of Solar system-sized highly 
inclined disks as faint as lOmJy in the mid-infrared. 

In this article we present a 5.2-37.2 //m spectrum of the edge-on disk CRBR 2422.8- 
3423 obtained with the Spitzer Space Telescope (Werner et al. 2004). CRBR 2422.8-3423 
(a =16''27'"28!8, 6 =-24°41'03"[J2000]) is a low luminosity source (~ 1 L©) located in the 
direction of the dense core Rho Oph-F in the Ophiuchus molecular cloud complex (Motte 
et al. 1998). The core contains a number of other young stars as well as two prestellar 
condensations. The brightest submillimeter point source in the core is IRS 43, located 34" 
to the north-east of CRBR 2422.8-3423. Near-infrared imaging has shown CRBR 2422.8- 
3423 to be a compact (1.5") bipolar nebula, characteristic of a disk viewed close to edge-on. 
The south-eastern lobe is a factor of 11 times fainter than the north-western lobe in the 
X-band (Brandner et al. 2000), showing that the disk has an inclination of ~ 70° (Thi et al. 
2002). The disk is unusually bright in the mid-infrared for an edge-on disk, making it an 
excellent target for infrared spectroscopy of ices. 

Previous studies of ices in disks have focused on the absorption band of solid CO at 
4.67 /xm observed toward isolated disks. These include Ehas 18 (Shuping et al. 2001), L 1489 
IRS (Boogert et al. 2002a) and CRBR 2422.8-3423 (Thi et al. 2002). Of these sources, the 
most inclined disk is probably that of CRBR 2422.8-3423, as evidenced in the near-infrared 
imaging of Brandner et al. (2000). Spitzer spectra have also been presented toward L1489 
IRS and DG Tau B (Watson et al. 2004). Common for all sources is the very red (increasing 
flux with wavelength) near-infrared colors indicating the presence of significant amounts of 
foreground or envelope material; an isolated edge-on disk is expected to have near-infrared 
colors that are much less red due to the dominance of scattering. The CO ice profiles along 
the three lines of sight are quite different, the first two being broad, indicating that much of 
the CO is embedded in a water-rich ice, while CRBR 2422.8-3423 is dominated by very deep 
absorption from largely pure CO ice (Pontoppidan et al. 2003). So far, only circumstantial 
evidence has been presented arguing for the location of the CO ice in the disks rather than in 
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envelope/foreground material. A key point in the observational study of ices in circumstellar 
disks is therefore to separate the contribution to the total ice bands of the disk from the 
contribution of any cold foreground or remnant envelope material. A detailed study of the 
location of ices toward the young star Elias 29 (located only a few arcminutes from CRBR 
2422.8-3423) found that most of the absorbing material along this particular line of sight 
is located in foreground clouds (Boogert et al. 2002b). This study provides an interesting 
reference to the case of CRBR 2422.8-3423. First, Ehas 29 is viewed face-on in contrast to 
CRBR 2422.8-3423. Second, it is has a bolometric luminosity of 36 Lq, a factor of 20-30 
higher than that of CRBR 2422.8-3423. Boogert et al. (2002b) find from the shallow depth 
of the CO ice band toward Elias 29 that the average depletion of gas-phase CO along this 
line of sight is only a few percent. 

This article will show that detailed radiative transfer modeling of high resolution near- 
infrared and/or millimeter imaging combined with high-quality mid-infrared spectroscopy 
can yield significant new information on the location and environment of ices. We use 
CRBR 2422.8-3423 as a case study, but emphasize that the method can easily be applied to 
other ~edge-on disks, using sensitive Spitzcr-IRS spectroscopy. This article is organized as 
follows: In §2, the observational data of CRBR 2422.8-3423 and the surrounding cloud are 
presented. §3 discusses the solid state absorption features detected in the infrared spectra. 
The radiative transfer model used for CRBR 2422.8-3423 is presented in §4. §5 discusses 
the distribution of ices in the CRBR 2422.8-3423 disk and §6 summarizes predictions of the 
model regarding observations of ices in other edge-on disks. 



2. Observations 

Observations of CRBR 2422.8-3423 at wavelengths from 1.2 /xm to 3 mm have been 
assembled using both new and archival data. JHKg imaging by the Infrared Spectrometer 
and Array Camera (ISAAC) mounted on UTl-Antu of the Very Large Telescope (VLT) 
has been extracted from the public ESO archive^. The ISAAC images were first presented 
by Brandner et al. (2000). L- and M-band spectroscopy has been obtained using the Long 
Wavelength spectroscopic mode of ISAAC; the M-band spectrum was discussed in (Thi et al. 
2002) ^, but we present a new observation of the ^'^CO ice band centered at 4.78 fim. Spitzcr 
5.2-37.2 iim spectroscopy was taken as part of the 'From Molecular Cores to Planet-forming 



^The VLT-ISAAC images were obtained at Paranal, Chile as part of the observing program 63.I-0691(A). 

^The VLT-ISAAC spectroscopy was obtained at Paranal, Chile as part of the observing programs 164.1- 
0605(A) and 71.C-0338(A) 
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Disks' (c2d) legacy program (Evans et al. 2003). A 350 /^m image was obtained with the 
Submilhmeter High Angular Resolution Camera II (SHARC-II) at the Caltech Submillimeter 
Observatory (CSO), and a SCUBA 850 /um map obtained with the James Clerk Maxwell 
Telescope (JCMT) was extracted from the 'COMPLETE survey of star-forming regions'^. 
Finally, a 5-16 /im spectral image obtained with the set of Circular Variable Filters on the 
Infrared Space Observatory (ISOCAM-CVF) was extracted from the public ISO archive^. 
The spectral image contains a number of sources in the Oph-F region and provides a useful 
comparison for the ice bands observed toward CRBR 2422.8-3423. 



2.1. Mid-infrcired spectroscopy 

5.2-37.2 fira spectroscopy was obtained with the Infrared Spectrograph (IRS) aboard 
the Spitzer Space Telescope on March 25, 2004 using the short-low (SL) module in the 
5.2-14.5 /im range with a spectral resolving power oi R — A/AA ~ 100, the short-high (SH) 
module at 9.9-19.6 //m and the long-high (LH) module at 18.7-37.2 //m with a resolving power 
of ~ 600. The archival AOR key is 0009346048 for PROGID 172. Total exposure times were 
28 s, 244 s and 484 s in the SL, SH and LH modules, respectively. The spectra were extracted 
from the pipeline (version S9.5.0) images using our c2d reduction software. The short-low 
spectra were extracted from pipeline products that were not corrected for stray light. The 
background in the SL modules was determined by fitting a Gaussian with an additive bias to 
the spectrum at each point along the dispersion direction. Special care was taken to match 
the orders of the SH and LH echelle spectra and remove the many bad pixels of the LH 
module. The SH module was scaled to the SL module in their overlap region, and the LH 
module was then scaled to the SH module. The final spectrum was smoothed with a 2-pixel 
box-car function. It should be stressed, however, that the present level of data reduction is 
likely to be improved significantly in the near future, especially in terms of the shape of solid 
state features. 

The L-band spectrum was obtained with ISAAC on the night of May 2, 2002 using the 
low-resolution mode and the 0'.'6 sht, yielding a resolving power oi R — 600. The R = 5 000 
M-band setting centered on the 4.67 ^um band of solid CO is dicussed in Thi et al. (2002). 
An additional M-band setting was obtained on May 8, 2003 of the 4.78 /xm band of solid 



^See http: / /cfa-www.harvard.edu/COMPLETE/ 

*In part based on observations with ISO, an ESA project with instruments funded by ESA Member States 
(especially the PI countries: France, Germany, the Netherlands and the United Kingdom.) and with the 
participation of ISAS and NASA. 



- 7- 



^^CO with the 0f!3 sht yielding R^IO 000. The ISAAC spectra were reduced following the 
procedure described in Pontoppidan et al. (2003). 

The ISOCAM-CVF image (AOR key 29601813) of the Oph-F core was reduced using 
the CAM Interactive Analysis (CIA) package ver. DECOl ^. A spectrum with it! = 35 was 
extracted from the 3x3pixels=18" centered on IRS 43, the nearest bright source to CRBR 
2422.8-3423. The CVF data were first presented by Alexander et al. (2003). 



2.2. (Sub)millimeter imaging 

The SHARC-II 350 //m and SCUBA 850 //m maps are shown in Fig. 1, overplotted on 
the VLT-ISAAC X-band image of the region. It is evident from the thermal dust emission 
that CRBR 2422.8-3423 is situated close to several embedded sources, at least in projection. 
IRS 43 is the brightest submillimeter source in the map. although the emission peak is 
shifted to the west by almost 10" in both maps. A significant extended source corresponding 
to the prestellar core Oph-F MM2 (Motte et al. 1998), is evident to the north-west of CRBR 
2422.8-3423. No obvious submilUmeter point source is associated with CRBR 2422.8-3423, 
although some emission is present slightly offset to the NW. It is a priori not clear whether 
the observed dust emission is due to emission from the edge-on disk, blended with extended 
emission or due to remnant envelope material surrounding CRBR 2422.8-3423. The offset 
emission from IRS 43 and CRBR 2422.8-3423 could indicate pointing errors. However, the 
offsets are similar in both the CSO and SCUBA maps. Additionally, IRS 44, IRS 46 and GY 
262 have 850 /xm point source emission, whose positions closely coincide with the infrared 
sources (Fig. 1). 

The disk is not detected at 3 mm with the Owens Valley millimeter interferometer with 
a 3(7 upper limit of 3mJy (8.5" x 4.0" synthesized beam). This upper limit to the amount of 
compact millimeter emission from the region is sufficient to indicate that the submillimeter 
emission is largely due to the core Oph-F MM2 together with possible extended remnant 
envelope material and not directly related to the disk. 



^The ISOCAM data presented in this article were analysed using "CIA", a joint development by the ESA 
Astrophysics division and the ISOCAM Consortium. The ISOCAM Consortium is led by the ISOCAM PI. 
C. Cesarsky. 
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Fig. 1.— Left panel: SHARC II 350 /im map of the CRBR 2422.8-3423 region overlaid on 
the ISAAC i^s-band image. The box indicates the area observed at 350 /xm. Contours start 
at 1.0 Jybeam"^ and increase in steps of 2.0 Jybeam"^. Right panel: SCUBA 850 //m map 
from the COMPLETE survey. Contours start at 0.05Jybeam~^ and increase in steps of 
0.05 Jybeam"^. North is up and east is to the left. 
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3. Inventory of ices 

The full 2.85-37.2 spectrum toward CRBR 2422.8-3423 is shown in Fig. 2. There 
are two holes in the spectral coverage at 4. 15-4.55 /xm and 4.9-5.2 /im. The first interval 
contains the CO2 stretching mode, which probably has a depth comparable to that of the 
CO stretching band at 4.67 /xm toward this source. The second interval is mostly observable 
from the ground, but it was not included in the ISAAC observations and is not expected to 
contain any strong solid state features. 

The spectrum is dominated by absorption bands due to sihcates and ices along the 
fine of sight. The strongest bands are those due to CO, H2O and CO2, but several weaker 
bands are also detected. Specifically, the 6.85 /im band, which still lacks an unambiguous 
identification, is prominent. In the following, the mid-infrared spectral features are discussed 
in order of increasing wavelength. Table 1 summarizes the abundances of the different ice 
species. 
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Fig. 2.— Combined VLT-ISAAC and Spitzer-IRS spectrum of CRBR 2422.8-3423. Detected 
solid state bands are marked. The water ice libration band at 11-13 //m is largely missing 
(see text). 



- 11 - 



3.1. X-band (2.85-4.15 /urn) 

The L-band spectrum shows a moderately deep 3.08 fim water ice band. The band is 
not saturated and therefore has a well-determined optical depth of 2.3 ± 0.3. The 3.54 fim 
band due to the C-H stretch of solid methanol (CH3OH) is not detected with an upper 
limit to the optical depth of 0.05, corresponding to an abundance relative to water ice of at 
most 10%. Similarly, the unidentified '3.47 //m feature', possibly due to NH3-H2O complexes 
(Dartois et al. 2002), is not detected. However, the hmits on the abundances are not very 
strict owing to the moderate water ice column density and the low flux levels of the L-band 
spectrum. Finally, there is no detection of hydrogen recombination lines such as Bra at 
4.07 /xm, possibly due to veiling by a strong near-infrared excess. 

3.2. M-band (4.55-4.90 /Lzm) 

The M-band spectrum shows an extremely deep CO ice band at 4.67 //m as well as 
narrow absorption lines from low- J rovibrational transitions of gaseous CO (Thi et al. 2002; 
Pontoppidan et al. 2003). In general, the shape of the main ^^CO ice band observed in dense 
clouds can be decomposed into three components with a fundamental set of centers and 
widths (Pontoppidan et al. 2003), such that all observed ^^CO ice bands can be well-fitted 
by only varying the relative intensities of the bands. The three-component decomposition is 
shown in Fig. 3. The CO ice band is dominated by the 'middle' component at 2139.9 cm~^, 
hkely due to mostly pure CO. The 'red' CO ice component at 2136.5 cm~^, probably due to 
CO in a water-rich environment, is detected with an optical depth of 10% of that of pure CO. 
The contribution of the 'blue' component centered on 2143.7 cm~^ is negligible toward CRBR 
2422.8-3423. This component is not unambiguously identified, but has been associated with 
CO mixed with CO2 (Boogert et al. 2002a) or with the longitudinal optical mode of crystalline 
CO, which appears when the background source is polarized (Pontoppidan et al. 2003). The 
dominance of the 2139.9 cm^^ component corresponds to about 80% of the CO molecules 
being in a mostly apolar environment. 

The ^'^CO band at 4.78 /xm is clearly detected with an optical depth of 0.17 ±0.03. This 
is only the third object for which ^^CO ice is found, and clearly testifies to the exceptionally 
large solid CO column density. In Fig. 4, the band is compared with laboratory spectroscopy 
from the Leiden database (Ehrenfreund et al. 1996). The -'^^CO band is useful for determining 
the molecular environment of the CO ice. This is due to the insensitivity of the ^^CO band to 
the grain shape effects that dominate that for ^^CO. Pure CO produces a narrow, Lorentzian 
^^CO band centered on 2092 cm~^, while relatively small contamination (at least 10-20%) 
from other molecules tends to broaden the band significantly. The presence of telluric CO 
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gas phase lines on either side of the ^^CO ice band, however, comphcates the profile analysis. 
Visual inspection of Fig. 4 indicates that the ^^CO band may be slightly broader than what 
is expected for pure CO. Assuming pure CO results in a ^^CO/^^CO ratio of 45 ± 15. This 
is somewhat lower than the value of 70 derived for other sources (Boogert ct al. 2002a; 
Pontoppidan et al. 2003). A broader ^"^CO ice band will result in a larger column density of 
^^CO and therefore an even lower isotopic ratio. 
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Fig. 3.— The ^^CO ice band observed toward CRBR 2422.8-3423 with VLT-ISAAC. The 
curves indicate the three-component decomposition of Pontoppidan et aL (2003). 
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Fig. 4.— The ^^CO ice band observed toward CRBR 2422.8-3423 with VLT-ISAAC. The 
dotted curve is pure CO at 10 K, the dashed and dashed-dotted curves are CO:CO2=l:0.7 and 
CO:CO2=l:0.21, respectively, both at 30 K. The difference in temperature of the laboratory 
spectra is not significant. Regions strongly affected by residuals from telluric lines have been 
removed from the spectrum. 
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3.3. The 5—12 /zm region 

The short-low region of the Spitzer spectrum, covering 5.2-14.5 ;um, contains clear de- 
tections of the enigmatic 6.0/6.85 fim solid-state band complex as well as the 9.7 fim silicate 
band. The silicate absorption band has an optical depth of 1.4-1.8, depending on how the 
continuum is defined, giving a small TH2o/''Siiicate ratio of ~1.4. This is often seen for embed- 
ded sources and may indicate that the silicate band has been partially filled in by emission 
(Alexander et al. 2003). 

Fig. 5 shows the 5-8 /xm region of the CRBR 2422.8-3423 Spitzer spectrum on an optical 
depth scale. Laboratory spectra (e.g. Schutte et al. 1996, 1999) as well the ISOCAM-CVF 
spectrum of IRS 43 and the ISO Short Wavelength Spectrometer (SWS) spectrum of the 
massive young star Mon R2 IRS 3 from Keane et al. (2001) are used for comparison. The 
continuum is a 2nd order polynomial chosen to pass through the 5.3-5.5 and the 14.8-14.9 
regions. There is some uncertainty involved in choosing a continuum in this way, but the 
shapes of the bands are not affected, even if their total optical depths may be uncertain by 
up to ~ 20%. The spectrum shows a prominent 6.0 fim feature, which is partly due to water 
ice. A laboratory spectrum of pure water ice deposited at 10 K and scaled to the 3.08 //.m 
band is compared to the observed spectrum. It is seen that the red wing of the 6.0 /xm 
water band can be fitted, while the main laboratory band can account for only 60% of the 
total observed band depth. This is very common and most likely indicates that other ice 
species contribute significantly to the 6.0 /xm band (Keane et al. 2001), although radiative 
transfer effects can also cause a similar effect (see §6). A significant blue shoulder to the 
6.0 fim band is detected at 5.85 /xm. This shoulder has been found in other high mass (Keane 
et al. 2001) and low-mass (Pontoppidan et al. 2004) sources. Plausible candidates include 
both formaldehyde (H2CO) and formic acid (HCOOH). In this case, HCOOH provides the 
best fit, both because pure H2CO is centered at 5.82 /xm, which is too blue for the observed 
spectrum, and because the 3.47 //m band of H2CO is not detected. A secure identification is 
not possible at the present, however, thanks to signal-to-noise limitations and level of data 
reduction. The HCOOH band shown in Fig. 5 corresponds to an abundance relative to 
water ice of 2.4%. 

A narrow absorption feature is seen at 6.35 /im. A sharp feature at this wavelength has 
not been observed for other sources. The band is possibly an instrumental artifact, although 
no obvious problems are seen in the pipeline spectroscopic images. Furthermore, the Spitzer 
Science Center pipeline extraction software ver. 9.5.0 also produces an absorption band at 
this wavelength comparable to that obtained with our own extraction. We therefore regard 
the detection of this band as tentative because it has not previously been reported toward any 
other source. A possible identification can be the formate ion (HCOO~), which has a strong 
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band centered around 6.33 //m. The presence of HCOO" is possible, given the presence of 
HCOOH, since it can be produced by acid-base reactions between, for instance, HCOOH 
and NH3. Only 0.4% of HCOO~ relative to water ice is required to produce a similar 
band. Another possible explanation is that an emission band due to polycyclic aromatic 
hydrocarbons (PAHs) is present at 6.2 fxm, thereby creating the illusion of an absorption 
band at 6.35 /xm. The spectral resolution is too low to distinguish between these possibilities. 
Finally, solid O2 has a very weak band at 6.45 //m. However, the depth of the observed band 
requires an extremely high column density of solid O2 of > lO^^cm"^ corresponding to an 
abundance of 2.5 x 10~^ for the largest band strength measured of 1 x 10~^^ cmmolec"^ 
(Ehrenfreund et al. 1992). The band strength of the O2 band has been measured to be 1-2 
orders of magnitude smaller, depending on ice composition (Vandenbussche et al. 1999). 
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Fig. 5. — The 5-8 /iin region observed toward CRBR 2422.8-3423. Upper panel: The observed 
spectrum is compared to a laboratory spectrum of water ice at 10 K scaled to an optical depth 
of 2.3 for the 3.1 /im band (dotted curve). A laboratory spectrum of pure HCOOH at 10 K 
has been added to account for the blue wing of the band (dashed curve). The solid curve 
is the sum of the two laboratory spectra. Lower panel: The spectrum is compared to the 
ISOCAM-CVF spectrum of the nearby IRS 43 (dotted curve) and the ISO-SWS spectrum of 
Mon R2 IRS 3 from Keane et al. (2001) (solid grey curve). Note that the red shoulder of the 
6.85 yum band is absent toward IRS 43 but can be well fitted by the Mon R2 IRS 3 spectrum. 
Three vertical dotted lines are added to guide the eye at 6.0, 6.85 /^m and 7.7 fim. The IRS 
43 spectrum has not been scaled to match the CRBR 2422.8-3423 spectrum, but happens to 
have almost the same optical depth. The dot-dashed curve is the radiative transfer model 
described in 54. 
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CRBR 2422.8-3423 also shows a prominent 6.85 //m band. The carrier of this band 
remains unknown, although a significant number of candidates have been tested. Substan- 
tial arguments have been given in favour of NH^ (Schutte & Khanna 2003, and references 
therein). A significant concern for a charged carrier of such a strong band is the lack of 
observed counter-ions, but Schutte & Khanna (2003) argue that the counter-ions may be 
distributed among a range of carriers, rendering any single ionic species unobservable. Keane 
et al. (2001) noted the presence of two components to the 6.85 //m band through an empirical 
decomposition of the band for a sample of high-mass stars. In Fig. 5 the 5-8 /im spectrum of 
CRBR 2422.8-3423 is compared to laboratory spectra as well as the ISOCAM-CVF spectrum 
of IRS 43 and the ISO-SWS spectrum of the massive young star Mon R2 IRS 3 from Keane 
ct al. (2001). Mon R2 IRS 3 is included here because it has the most extreme red-shifted 
6.85 yum band observed to date. The 6.0/6.85 //m band complexes toward CRBR 2422.8- 
3423 and IRS 43 have, by chance, nearly identical optical depths, while the spectrum of Mon 
R2 IRS 3 has been scaled by a factor of 0.7. The main difference is the complete lack of 
the red shoulder of the 6.85 /xm band toward IRS 43 as compared to CRBR 2422.8-3423. 
This difference may be of considerable significance. The Mon R2 IRS 3 spectrum is seen to 
provide a very good fit to the red part of the CRBR 2422.8-3423 6.85 //m band, which is 
entirely missing from the IRS 43 spectrum. A shift of the 6.85 /im component to 7.0 /im was 
suggested by Keane ct al. (2001) to be an effect of thermal processing, meaning that the line 
of sight toward CRBR 2422.8-3423 contains warm ice in addition to the cold ice probed by 
the CO ice band. Schutte & Khanna (2003) have shown that NH^ is shifted to the position 
of the reddest observed band toward Mon R2 IRS 3 at temperatures in excess of 200 K. The 
corresponding temperature for interstellar ice is possibly somewhat lower. The location of 
material at high temperatures along the hne of sight toward CRBR 2422.8-3423 is discussed 
in terms of the radiative transfer modeling in §5. The abundance of NH4 required to produce 
a 6.85 iim feature like that observed toward CRBR 2422.8-3423 is 10% relative to water ice, 
assuming a band strength of 4.0 x 10~^^cmmolec~^ (Schutte & Khanna 2003). 

The 7.7 /im band of solid CH4 is tentatively detected, but unresolved at the Spitzer-IRS 
resolution. The silicate absorption band at 9.7 /xm is largely featureless (see Fig. 2) and does 
not show strong evidence for the 9.35 iim 'umbrella' mode due to solid NH3 or the 9.75 iim 
CO stretching mode of solid CH3OH. 



3.4. CO2 15.2 /xm bending mode 



The short-high spectrum shows a strong band at 15.2 //m due to the bending mode of 
CO2 ice. The CO2 bending mode observed toward CRBR 2422.8-3423 is shown in Fig. 6 
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where it is compared to the laboratory spectroscopy of Ehrenfreund et al. (1996) as well as 
the CO2 bending mode observed toward IRS 43. The depth of the band corresponds to an 
abundance of CO2 ice relative to water ice of iV(C02)/iV(H20) = 0.32 (see Table 1). This is 
significantly higher than the typical value of 0.2 observed for high-mass stars (Gerakines et al. 
1999), but is similar to abundances recently measured toward some low-mass stars (Boogert 
et al. 2004). Laboratory experiments have shown that the shape of the CO2 bending mode is 
a sensitive tracer of ice composition and thermal history. In particular, substructure bands 
can be used to examine whether solid methanol is mixed with the CO2 and the extent of 
ice heating (Ehrenfreund et al. 1999). The location of the CO2 bending mode within the 
Spitzer high spectral resolution coverage enables a detailed study of the band shape. 
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Table 1: Average abundances of ices observed toward CRBR 2422.8-3423 





Column density 


Abun. rel. to H2O ice 


Abun. rel. to H2 gas 


band strength 


H20 


3.6 


1 


9 X 10-5 


2 X 10-16 a 


C02 


1.16 


0.32 


2.9 X 10-5 


1.1 X 10-1^" 


12 CO pure 


2.8 


0.78 


7 X 10-5 


1.1 X 10-1^" 


12 CO in water 


0.68 


0.19 


1.7 X 10-5 


1.1 X 10-1^" 




0.36 


0.10 


9 X 10-*5 


4.0 X 10-^"^" 


"CO 


0.06 


0.017 


1.5 X 10-6 


1.1 X 10-1^"* 




0.13 


0.035 


3 X 10-6 


8.0 X 10-1^^ 


CH3OH 


< 0.36 


< 0.10 


< 9 X 10-6 


5.3 X 10-1^^ 


HCOOH'' 


0.09 


0.02 


2 X 10-6 


6.7 X 10-1^^ 



" Gerakines et al. (1995) 



Kerkhof et al. (1999) 
^ Schutte et al. (1999) 

Suggested identification (see text). 
^ Schutte & Khanna (2003) 
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Pure, partially crystalline CO2 at 10 K is split into two narrow bands at 660 cm~^ 
(15.15 /xm) and 655 cm~^ (15.25 /im). There is some evidence from laboratory experiments 
that pure, fully amorphous CO2, which does not show this characteristic splitting of the 
bending mode, may be formed under certain conditions (Falk 1987). Toward CRBR 2422.8- 
3423, a strong, relatively narrow component is clearly present around 660 cm~^ (15.15;um). 
Inspection of Fig. 6 shows that a narrow band at 655 cm~^ is absent, ruling out significant 
amounts of the pure, partially crystaUine CO2 toward CRBR 2422.8-3423. The absence of the 
655 cm~^ band indicates that the structure of the CO2 ice is different from that generally ob- 
served toward high-mass stars (Gerakines et al. 1999). In addition to the narrow component, 
the observed CO2 band has a significant red shoulder at 647 cm~^ (15.45 //m), which can pos- 
sibly be assigned to a CO2-CH3OH complex (Gerakines et al. 1999). This shoulder is most 
prominent for cold ices (< 50 K) and when CO2, H2O and CH3OH are present in the mix- 
ture in similar concentrations. A very good fit is obtained with a H2O:CH3OH:CO2=l:0.6:l 
mixture. The abundance of CH3OH required by the fit is 9% relative to water ice, which is 
close to, but consistent with the observed upper limit of methanol ice. 

A better theoretical and experimental understanding of the conditions under which 
different phases of solid CO2 are formed will allow us to explore the structure of CO2 ice 
around low-mass protostars further. For example, an alternative way to create a single 
narrow absorption band of GO2 is to use a mixture with another ice species. Since a mixture 
with H2O cannot produce the narrow blue component, the only other candidate known to 
be abundant in this source is GO. Indeed, a binary mixture with GO is found to yield an 
excellent fit to the rest of the band, and the percentage of the total column of 'pure' GO 
required to be mixed with GO2 is 20%. 

The presence of a binary G0:G02 component is quahtatively consistent with the ^^GO 
band. Quantitatively, the best fit to the GO2 bending mode is a mixture with GO and GO2 
in approximately equal concentration, while the ^^GO band is best fitted by a mixture in 
which GO2 is much more dilute. Since the total ratio of CO to CO2 along the line of sight 
is ~ 3, this fits well with a scenario in which a small part (~ 20%) of the CO is mixed with 
the CO2, while the rest is almost pure CO. A binary C0:C02 mixture should produce a 
blue-shifted and broadened ^^CO ice band, and there is no indication of this in the observed 
M-band spectrum (see §3.2). Since the GO ice band toward CRBR 2422.8-3423 is saturated, 
the significance of this is not clear. 

In summary, the line of sight toward CRBR 2422.8-3423 shows significant evidence for 
some interaction between part of the CO2 and the CO ice components, while most of the 
two species remain separate. Alternatively, the narrow component of the GO2 bending mode 
may be due to pure, but amorphous GO2. One may speculate that the G0:G02 ice mixture 
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or amorphous CO2 ice resides in the disk material, mostly on the basis of circumstantial 
evidence. First, this component seems to be much less prominent toward IRS 43, although 
an observation of this line of sight at higher spectral resolution is required to confirm this. 
Second, a broadened ^'^CO band has not previously been observed in lines of sight toward 
embedded sources. However, considering that CRBR 2422.8-3423 is only the third source 
with a ^^CO detection, this can hardly be considered conclusive. 
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Fig. 6. — The 15.2 jjm CO2 ice bending mode observed toward CRBR 2422.8-3422 on an 
optical depth scale. Upper panel: The curves show a sum of laboratory spectra of a water 
and methanol-rich mixture (H2O:CH3OH:CO2=l:0.6:l at 30 K) and a CO rich mixture. 
The solid curve shows a fit using CO:CO2=l:0.7 at 10 K for the CO-rich mixture, while 
the dashed curve uses CO:CO2=l:0.2 at 10 K. A continuous distribution of ellipsoid (CDE) 
grain model has been used to convert the optical constants to absorption coefficients. Lower 
panel: The dotted curve is the CO2 band observed by ISOCAM-CVF toward IRS 43. The 
dot-dashed curve is the band extracted from our radiative transfer model (see §4). The 
spectral resolution of the ISOCAM spectrum is 0.4 /im at the location of the CO2 bending 
mode. 
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4. 2D continuum radiative transfer 



In order to better understand the location of ices in the disk and the foreground or 
envelope of CRBR 2422.8-3423, a 3D Monte Carlo radiative transfer model was constructed 
for an axisymmetric density structure. The aim was to fit the dust density structure simulta- 
neously with the high resolution spectral energy distribution and the near-infrared imaging. 
We use the RADMC code (see DuUemond & Dominik 2004) to calculate the temperature 
structure and scattering source function for a given density distribution. Spectral energy 
distributions and images are then calculated using the ray-tracing capabilities of the code 
RADICAL (Dullemond & Turolla 2000). The Monte Carlo radiative transfer is effectively 
performed in full 3D, but the density structure is required to be axisymmetric. The scat- 
tering is assumed to be isotropic, although a full treatment of scattering is being developed. 
The effects from anisotropic scattering are relatively small for wavelengths longer than 2 /^m 
and grains smaller than a few /im. However, this does not necessarily rule out that the 
near-infrared images are somewhat affected by anisotropic scattering. 



4.1. Dust model 



Because it is necessary to adjust the abundances of ice species independently, a dust 
model must be constructed for the radiative transfer calculation to provide opacities that 
can reproduce the observed spectral features. Furthermore, to adequately fit the relatively 
narrow ice bands, a higher spectral resolution is required for the opacities than is available in 
published tables. The adopted dust model assumes spherical silicate grains with inclusions 
of carbonaceous material, surrounded by a spherical ice mantle. In the diffuse medium, 
the silicate and carbonaceous grains are likely in separate dust populations (Draine 2003). 
However, when grains coagulate in dense molecular clouds, these populations may be mixed. 
Thus for convenience, a single grain population is adopted here. The silicate-carbon core 
optical constants are calculated using Maxwell-Garnett effective medium theory (Bohren & 
Huffman 1983): 



1 + 



3/( Cine ^mat)/ (^inc ~l~ '^^uiat) 



1) 



1 ~ /(^inc ~ ^mat)/ (^inc + '^^ma.t) 

where emat ^-nd einc are the dielectric functions of the matrix (silicates) and spherical inclusions 
(carbon), respectively. / is the volume fraction of the inclusions. The carbon component 
is chosen as the inclusion, since it has the smallest volume fraction. The opacities are then 
obtained using the coated sphere Mie code from Bohren & Huffman (1983). A single grain 
size of 0.5 /im is used, since the shape of the opacity law is mostly affected by smaller grains 
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at the shorter wavelengths (< 2 /im), where httle information is available due to the heavy 
extinction of the source. Ices are assumed to evaporate at 90 K, except CO, which evaporates 
at 20 K (Sandford & Allamandola 1993). Above 90 K, only the silicate/carbon cores are left. 
The grains are assumed to be compact. 'Fluffy', coagulated grains have increased opacity at 
all wavelengths (Ossenkopf & Henning 1994); consequently, we will overestimate the required 
dust mass by assuming compact grains. 

Optical constants for the ice are taken from the Leiden database of ices (Ehrenfreund 
et al. 1997), those for the oxygen-rich silicates from Ossenkopf & Henning (1994), while the 
carbonaceous constants arc from Jagcr ct al. (1998) for carbon clusters formed at 800 K. The 
choice of the carbonaceous component is the most uncertain, since few spectral diagnostics 
are available in dense clouds and embedded objects to guide the choice. The carbonaceous 
material created in the laboratory experiences a significant change in properties for pyroliza- 
tion temperatures between 600 and 800 K, reflecting a larger contribution from graphitic 
material at the higher temperature. Here, the optical constants for carbon dust created at 
800 K arc chosen because this best reproduces the observed strength of the 9.7/im feature. 
The 600 K sample results in opacities that cause the 9.7/im silicate feature to be over- 
estimated in the model fit to the observed spectrum, while creating very high far-infrared 
to millimeter opacities. The volume fraction of carbonaceous inclusions used is 0.15. This 
number is dependent on the chosen set of optical constants for the carbon dust. A lower 
pyrolization temperature creates a lower carbon dust opacity in the mid-infrared, causing 
the volume fraction of carbon dust required to fit the silicate feature to increase by up to 
a factor of two. The significance of the carbon volume fraction is therefore hmited. For 
comparison, the dust model of Weingartner & Draine (2001) uses a carbon volume fraction 
of ~ 0.3. 

It is not possible to create a database of ice optical constants that covers even a small 
fraction of the parameter space of ice abundances. As a consequence, ice abundances in the 
mantle are varied using Eq. 1. As shown in §3, much of the CO2 ice is embedded in water 
ice. Therefore, as a starting point an ice mixture of H2O:CO2:CO=l:0.2:0.03 at 10 K is used. 
The CO2 and CO abundance is increased by adding inclusions of CO:CO2=l:0.7. Beyond 
20 iim, a revised version of the ice opacity of Warren (1984) for crystalline water ice is used, 
since no optical constants are available for amorphous water ice at these wavelengths. The 
C02:C0 inclusions can be replaced by pure CO to test the location of this component in 
the disk/envelope system (see §5). The water ice abundance for the ice-coated dust grains 
is 9 X 10~^ relative to H2, assuming an ice density of 0.6gcm~^ and a gas-to-dust mass ratio 
of 100. 

Fig. 7 compares the dust opacities (not including the gas mass) to those by Ossenkopf 
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& Henning (1994) and Weingartner Sz Draine (2001). The calculated opacity matches fairly 
well that of Ossenkopf & Henning (1994) for compact bare grains (the first column of their 
Table 1), whose grain coagulation model produces opacities that are slightly larger at all 
wavelengths. The opacities show a larger variation once ice mantles are added. In particular, 
our higher resolution opacities include the narrow ice bands from CO and CO2 at 4.27 and 
4.67/xm. Additionally, we also find a significant contribution from the 11-13 /^.m water ice 
libration band in our opacity. Note also the presence of the 45 /im band from crystalline water 
ice due to the use of optical constants for higher temperature water ice at long wavelengths. 
The presence of this band does not affect the continuum modeling results. 
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Fig. 7. — Examples of dust mass absorption (top) and scattering (bottom) efficiencies used as 
input for the radiative transfer model compared to opacities from the literature. Left panels: 
Mass absorption and scattering coefficients for dust grains without ice. The fairly high 
absorption coefficients at long wavelengths are largely due to the choice of carbonaceous 
optical constants. Right panels: Dust mass absorption and scattering coefficients for ice- 
coated grains with a water ice abundance of 9 x 10^^ relative to H2, assuming a gas-to-dust 
ratio of 100 and an ice density of 0.6gcm~^ (porous amorphous water ice). The opacities are 
relative to the mass of the refractory dust mass in order to preserve the dust-to-gas mass in a 
model when also ice freezes out. The opacities are compared to those of Ossenkopf & Henning 
(1994) for compact grains with and without a thin ice mantle and those of Weingartner & 
Draine (2001) appropriate for diffuse medium dust. 
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4.2. 



Disk structure 



The aim of the radiative transfer model is to empirically fit the dust structure of the 
disk, rather than to construct a self-consistent dynamical model. Basically, this means that 
a 'generic' disk structure is used in which the flaring angle is a free parameter, independent 
of the central object or the dust model. The vertical density structure is assumed to be 
Gaussian. A Gaussian density profile is appropriate for a vertically isothermal disk. Note, 
however, that the model disk is not isothermal and no assumption is made concerning the 
temperature, which is calculated independently using the Monte Carlo code. Non-isothermal 
models have vertical density profiles that are similar to a Gaussian, only in the uppermost 
layers are noticable deviations seen (Dullemond et al. 2002). Effects such as grain settling 
will produce a vertical dust density profile which peaks more strongly toward the disk mid- 
plane, but here we are most interested in modeling the structure of the disk where the 
optical depth at mid-infrared wavelengths is moderate, away from the densest parts of the 
mid-plane. A dust density profile which is strongly peaked toward the mid-plane (e.g. due 
to dust sedimentation) is not probed by the mid-infrared spectrum as long as the wing of 
the profile resembles that of a Gaussian with a larger scale height, so that the photosphere 
of the disk lies well above the midplane density concentration. We therefore do not include 
a midplane dust layer into our model. The disk density is thus modeled as a function of 
radius, R and height above the disk plane, Z: 



is the disk scale height. The flaring power of 2/7 is that determined for a self- irradiated 
passive disk (Chiang & Goldreich 1997). The parameters varied in the model are Edisk, 
determining the disk mass, and -ffdisk determining the disk opening angle measured at the 
outer edge of the disk, -Rdisk- The disk opening angle can be constrained with some degree of 
confidence by the near- infrared imaging. The surface density is assumed to vary with p = 1. 
The inner edge of the disk is determined by the dust sublimation temperature, here assumed 
to be 1700 K, and a lower bound to i?disk is constrained by the X-band image to 90 AU (see 
Fig. 9). The surface density is not sharply truncated at the outer edge, but follows a steep 
power law of R~^'^ outside of -Rdisk- 





where S(i?) = Edisk x {R/ Rdisk) ^ is the surface density and 



Hp{R)/R— {H disk/ Rdisk) X {R/ Rdisk)'^^'^ 



(3) 



The central object is assumed to be a typical pre-main sequence star with an effective 
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temperature of ~ 4500 K. The radius of 2 solar radii (corresponding to a luminosity of 1.4 Lq) 

is then constrained by the SED. These parameters correspond to a pre-main sequence star 
of age 3 X 10^ years, using the evolutionary tracks of Siess et al. (2000). The model output 
is not very sensitive to the exact choice of temperature of the central star, only to the 
total luminosity. Specifically, the star could be cooler and therefore younger. An effective 
temperature of 4000 K, for example, corresponds to an age of 10^ years at the luminosity of 
CRBR 2422.8-3423. 

Foreground or envelope material is added via a spherical shell of constant density. This 
is necessary to produce the observed high extinction toward the source. The temperature 
and any deviation from a constant density of the foreground material is not constrained due 
to the lack of far-infrared photometry. The location of the foreground material is discussed 
in the next section. 

Finally, a puffed-up inner rim is added. The rationale behind a puffed-up inner rim 
is that the star over-heats an inner vertical wall of dust located at the dust sublimation 
temperature, causing it to expand in the vertical direction (Natta et al. 2001). This is 
modeled as a puffed-up scale height at the inner edge of the disk and a drop-off to a maximum 
radius where the flaring disk takes over. Specifically, the disk opening angle at the inner 
radius, i?rimj is increased to i^rim- The presence of a puffed-up inner rim makes it possible 
to reproduce the high mid-infrared fluxes of the source, and at the same time sufficiently 
extinct the central star. It is thought that a puffed-up inner rim is most prominent for 
disks surrounding intermediate mass stars (Natta et al. 2001; Dullemond et al. 2001), but 
it has also recently been found to be important for T Tauri stars (MuzeroUe et al. 2003). 
The observed SED of CRBR 2422.8-3423 clearly requires signiflcant excess infrared emission 
around 5 /im, which is hard to explain without introducing an inner rim. Also, at the 
inclination of CRBR 2422.8-3423, direct emission from the central star tends to dominate 
over scattered light from the disk surface if no inner rim is used, contrary to what is observed 
in the near-infrared images. 

The density distribution is sampled on a polar {R,Q) grid with a logarithmic spacing 
in the i?-direction to ensure that all structures are properly resolved. Likewise, the spectral 
energy distributions are calculated using rays with impact parameters that are sampled on 
a logarithmic grid centered on the star. At each wavelength, the flux is determined by 
integrating over a predefined 'observational' aperture. The images are calculated as follows: 
Pixels are arranged on a regular rectangular matrix and ray-tracing along the ray belonging 
to each pixel is performed. However, some extra care has to be taken with this procedure. 
Because the pixels are equi-distant, small scale structure in the inner parts of the disk (such 
as the star itself) may not be resolved. This will cause an underestimate of the flux in the 
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inner pixel of the grid. To ensure that the image reproduces the flux at all scales, including 
scales that are much smaller than a pixel, the following approach was taken: For each final 
model image, a set of three images was calculated on three different scales, all centered on the 
central star. The sampling chosen for the final image was 200 pixels for 1000 AU. Therefore, 
two images with 200 pixel sampling were calculated for the inner 10 AU and 0.1 AU. The 
smallest image was then resampled to a 2 x 2 pixel grid, which was used to replace the inner 
pixels in the 10 AU image. This method was then repeated for the 1000 AU image. The 
smallest image fully resolves the central star. 



4.3. Best model fit 
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Fig. 8. — Structure of the best-fitting model disk. The left panels show the entire disk, while 
the right panels show a zoom-in of the inner rim region. Upper panels: Dust temperature in 
Kelvin. The dashed curve indicates the (water) ice evaporation zone at 90 K. Lower panels: 
H2 density assuming a gas-to-dust ratio of 100 in units of cm~^. Contours are shown in steps 
of 1 dex (left) and 2 dex (right). 
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Fig. 9. — Observed and modeled composite J (blue), H (green), Kg (red) images of CRBR 
2422.8-3423. The color scale is logarithmic. The left side show the images of the disk as 
they appear including the foreground material, while the right side shows the images with 
the foreground material removed. Note that the central star and inner rim are seen as a 
reddened (by the disk material) point source in the center of the disk when the foreground 
is removed. The images are flux calibrated such that the observed and modeled colors can 
be directly compared. Lower panel: VLT-ISAAC color composite. Middle panels: JHKs 
model image corresponding to the SED in Fig. 12. The image has been convolved with 
a gaussian of FWHM = 0'.'4 to match the observed seeing. Upper panels: As above, but 
without convolving with the instrumental point spread function. 
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The 'best fit' density structure shown in Fig. 8 corresponds to a set of parameters 
which reproduce the observed constraints. Fig. 9 compares the synthetic JHKg composite 
images of the model disk to the observed image. Contours of the ISAAC images compared 
to the 'best fit' model are shown in Fig. 10, while contours showing alternative model 
images discussed in §4.4 are shown in Fig. 11. Figs. 12 and 13 compare the modelled 
and observed SED. The best-fitting parameters are summarized in Table 2. The table also 
includes estimates of the uncertainties of each model parameter due to degeneracies. 

The model can, in principle, accommodate an arbitrary number of dust opacities, which 
can be made to depend on the dust temperature. This is useful for simulating the evaporation 
of different ices at different temperatures. For instance, a layer of solid pure CO desorbs at 
20 K at densities of ~ lO^cm"^, while water ice desorbs at 90 K (Sandford & AUamandola 
1993). However, the desorption temperatures are very sensitive to the structure of the ice. 
CO may be trapped in water ice and desorb at much higher temperatures than pure CO 
(Collings et al. 2003), for example. Here we make a very simple model, assuming as a starting 
point that CO desorbs at 20 K, since it is known to be largely pure, while all other ice species 
desorb at 90K. 

The presence of several temperature-dependent opacities requires that the model is 
iterated to obtain a self-consistent temperature structure. This was done by first calculating 
a model with ice-covered grains throughout the disk. The output temperature structure 
was then used to determine where the ices evaporate and a new model was calculated with 
opacities uniquely defined at the relevant locations in the disk. This was repeated until the 
temperature became stable within 5%. It was found that convergence is usually reached in 
the first few iterations. 

The parameters were varied manually to obtain the best fit because computing time 
prohibits the calculation of a full grid. Any degeneracy in the best-fitting model parameters 
is therefore difficult to assess in detail. However, some parameters are clearly independent. 
As mentioned in §4.2, the size of the flared disk and the disk opening angle are constrained 
by the near-infrared images, as is the inclination angle. It is important to note that the 
disk may have a non-flared, self-shadowed component at large radii, which would not be 
seen in the scattered light in the near-infrared. Thus, the near-infrared images constrain 
the size of the flared disk, but provide only a strict lower limit to the size of the entire 
disk if a hypothetical non-fiared outer part is included. However, the current data cannot 
constrain a self-shadowed component. In particular, a fiat outer disk will not contribute to 
the observed ice bands since the line of sight to the infrared source does not pass through 
such a flat disk component. Therefore only the flared part of the disk is modelled. While 
not unambiguously determined, the total mass of the disk is constrained by the upper limit 
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observed at 3 mm in the OVRO interferometer to be < 0.005 Mq, using a dust mass opacity 
of K(3mm) = 0.08cm^g~^. This upper hmit also weakly constrains the possible mass of an 
outer self-shadowed disk component. 
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Fig. 10.— Left panels: Contour plots of the VLT-ISAAC JHKs images of CRBR 2422.8- 
3423. Right panels: Contour plots of the JHKs model images for the parameters given in 
Table 2. Note that the central star is visible in the model images. The lowest contour is at 
1.6 X 10~^^ergs~^ cm~^ Hz~^ sterad"^ and subsequent contours increase in steps of 0.4 dex. 
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Fig. 11. — Left panels: Contour plots of the JHKs model images for the 'standard' model, 
but with a small power-law envelope added to account for the excess nebulosity seen in 
the ISAAC imaging. Right panels: Contour plots of the JHKs model images for a disk 
with -Rdisk ~ 180 AU. The lowest contour is at 1.6 x lO^^^ergs"^ cm~^ Hz~^ sterad"^ and 
subsequent contours increase in steps of 0.4 dex. 
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A more realistic grain size distribution may significantly affect the mass estimate of the 
disk (Wood et al. 2002). However, the shape of the ice bands as well as the 9.7 fjm silicate 
band indicate that the grains contributing to the absorption bands are dominated by grain 
sizes of 0.1-1.0 ;um. This is due to the fact that scattering on grains with 27ra/Aicc ^ 1, 
where a is the size of the grain and Aicc is the wavelength of the ice band in question, 
will enhance the contribution to the band from the real part of the refractive index, thereby 
creating excess absorption on the red side of each band. This is typically seen for the 3.08 //m 
water ice band (see e.g. Leger et al. 1983; Smith et al. 1989; Dartois & d'Hendecourt 2001). 
For CRBR 2422.8-3423, the red 'scattering wings' of the ice bands are weak at wavelengths 
longer than 4 /im, indicating that the contribution from grains with sizes ^ 1 /xm is small. 
However, this may also simply indicate that the ice absorption is dominated by foreground 
material. 

The most important indirect constraint on the disk structure is the high level of near to 
mid-infrared emission. Given the inclination angle, a disk mass of 0.01 Mq or larger produces 
too much extinction along the line of sight to the inner disk to match the Spitzer spectrum. 
Even with the lower mass used in the presented model, it is found that the presence of an 
inner rim is essential to reproduce the shape and high flux level of the mid-infrared spectrum. 
Alternatively, unrealistically lowering the effective temperature of the central star to 2000 K 
produces a similar effect. An important point in this context concerns the visibility of the 
central star in the near-infrared images. At an inclination of ~ 90° x (1 — i?disk/-Rdisk), the 
inner disk and star become visible to the observer. This usually results in the star outshining 
the disk by orders of magnitude. This is clearly not the case for CRBR 2422.8-3423 and is 
the reason the source was first detected as an edge-on disk. However, the blocking of the star 
by the fiared disk also tends to strongly suppress the mid-infrared continuum. Therefore, in 
order to fit both the images as well as the Spitzer spectrum, only a very narrow range of 
surface density profiles, and consequently disk masses, arc possible, given an opacity. The 
puffed-up inner rim is vital for blocking the light from the central star sufficiently without 
suppressing the mid-infrared continuum, although even in our best fit model some of the 
stellar light can still be seen directly in excess of what is observed in the X-band; in the 
model the star is about 10 times too bright at 2.2 //m. However, in the J- and if-bands 
the model star is much fainter than the disk nebulosity in accordance with the observations. 
Furthermore, the bright mid-infrared continuum clearly indicates that the central parts of the 
disk dominate over scattered light at wavelengths longer than ~3 fim. Thus, the problem is 
confined to the A'-band image. This is most clear in Fig. 9 where the foreground component 
has been removed. Here a reddened (by the disk material) point source is visible in the center 
of the disk. The fact that the central star in the model cannot be completely suppressed 
at 2 /im while preserving the high mid-infrared continuum is a significant short-coming of 
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the model. Possible explanations include a different near-infrared opacity law or a different 

geometry of the inner disk, in particular the puffed-up inner rim. The brightness ratio of 
1:11 (Brandner et al. 2000) of the two scattering lobes of the disk was therefore fitted by 
comparing the modeled and observed i^-band images along a line offset by 50 AU along 
the major axis of the disk to avoid the light of the central star. The 'best-fitting' model 
brightness ratio is 1:12 along this offset line. 

The column density of foregroimd material toward the disk is well constrained by the 
near-infrared color to 2.2 x 10^^H2cm^^, assuming a gas-to-dust mass ratio of 100, corre- 
sponding to an extinction in the J-band of 8 mag. If this material is placed close to CRBR 
2422.8-3423, for instance in the form of a remnant envelope, a bright reffection nebulosity 
is created. Since no such extended reflection nebulosity is seen in the K-hand imaging, the 
foreground material must be located at a distance of at least 1000 AU from the source. This 
is consistent with the line of sight passing through the Oph-F core seen in the submillimeter 
imaging, rather than an envelope surrounding CRBR 2422.8-3423. Alternatively, the source 
may already have excavated a significant cavity in its envelope. A small amount of excess 
nebulosity compared to the 'best fit' model within 300 AU is seen in Fig. 10. The possible 
presence of a small tenuous envelope of CRBR 2422.8-3423 is discussed in §4.4. The lack of 
bright, extended scattering nebulosity in the near-infrared imaging puts an upper limit of a 
fewx 10^ for the density within 1000 AU of CRBR 2422.8-3423. The density of the material 
producing the near-infrared extinction was constrained by fitting the submillimeter emission 
within an aperture of 30" given the column density calculated using the near-infrared col- 
ors. The arbitrary method of fitting the foreground material using a spherically symmetric 
envelope introduces an uncertainty of a factor of 2 in the density. This is because the source 
may be located entirely behind the material producing the submillimeter emission and not 
in the middle, such as is assumed in the model. 
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Table 2: Best-fitting model parameters 





Actual model value 


Estimated range (see text) 




4500 K 


3500 5500 K 




2.0 i?o 


1.5-2.5 i?0 




0.8 M© 


0.2-2.0 M© 




IALq 


1.3-1.51.0 


Mdisk 


0.0015 Mq 


0.0005-0.005 M0 


-Rdisk 


90 AU 


80-100 AU (Only flared disk) 


-f^disk/ ^disk 


0.18 


0.16-0.20 


-^^rim/ Rrim 


0.12 


0.10-0.14 


R'm^env 


1000 AU 


> 1000 AU 


-^out,env 


6700 AU 


< 7000 AU 


Penv 


2.6 X lO^cm-3 


1.3-3.0 X lO^cm-3 


Incl. angle 


69.2° 


68-72° 


H2O ice abundance (env.+disk) 


9 X 10-^ (rel. to H2) 


8- 10 X 10-5 


CO ice abundance (envelope) 


1.3 X 10-4 (rel. to H2) 


Assuming no CO ice in the disk 
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Fig. 12.— Best fitting SED model of CRBR 2422.8-3423. The solid gray curve shows 
the model fiux through an aperture of 10", while the dashed line shows the fiux through 
an aperture of 30". The 350 //m, 850 iim and 1300 iim fiuxes are integrated though a 30" 
aperture centered on the near-infrared position of CRBR 2422.8-3423. The 1300 /xm fiux is 
taken from Motte et al. (1998). 
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Fig. 13. — A close-up comparison of the best-fitting model in observing apertures of 10" 
(grey curve) and 30" (dashed curve) with the observed infrared spectra (black curve). Note 
the mismatch on the red side of the 9.7/xm silicate feature due to the water libration mode 
in the model. 



4.4. Alternative models 



In this section, we explore the sensitivity of the model to simultaneous SED and imaging 
constraints. Specifically, we explore if the temperature of the absorbing material in the disk 
can be lower than in the 'best-fitting' model in order to facilitate the formation of CO ice. In 
this context, an important question is whether the size of the flared disk is well-constrained 
by the ISAAC image, since a larger disk may contain colder dust. To test this, a model was 
calculated by assuming i?disk = 180 AU, but flxing the remaining structural parameters to 
those in the 'standard' model given in Table 2. A best fit to the SED was then obtained 
by varying the luminosity of the central star and the inclination of the system. Because 
the mass of the disk is kept constant the inclination needs to be higher to produce the same 
extinction toward the star. At the same time, the fraction of warm dust in the disk decreases, 
requiring the luminosity to be increased to fit the mid-infrared emission. This uncertainty on 
the luminosity is primarily due to the lack of far-infrared photometry. The resulting image 
is seen in Fig. 11 and the SED is compared to the 'standard' model in Fig 14. It is seen 
that the SED of the large disk can still provide a good fit to the data by varying only the 
inclination and luminosity. However, the Ks image of the large disk clearly becomes larger 
than the observed image. We conclude that the size of the fiared disk is constrained by the 
near-infrared imaging. Modeling of an outer non-flared (self-shadowed) disk is postponed to 
a later paper. 

The ISAAC Ks image also clearly shows excess scattering emission close to the central 
disk (Fig. 10). This emission is not included in the 'standard' model. A small envelope may 
affect the modelled flaring properties of the disk by raising the scattering surface (Stapelfeldt 
et al. 2003). Here, we test the effect on the model image and SED of a small envelope 
producing excess emission. A small power-law envelope is added at radii larger than i?disk — 
90 AU to flt the ISAAC images. An envelope with a density power-law exponent of -1 and a 
gas density at -Rdisk of 10^ cm~^ provides a reasonable flt to the excess emission in the ISAAC 
Kg image. Note, however, that the shape of the density proflle is not strongly constrained. 
The influence of the small scattering envelope on the SED is seen in Fig. 14 to be minimal. 
Additionally, the column density of such a remnant envelope is so small that it has no 
influence on the ice absorption bands in the model. 
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Fig. 14.— Comparison of the SEDs for different models of CRBR 2422.8-3423. The dashed 
curve shows the 'standard' model fit presented in Fig. 12. The dash-dotted curve shows the 
SED if a small power-law envelope has been added to account for the excess nebulosity seen 
in the ISAAC images of CRBR 2422.8-3423. The sohd curve shows the SED of a model with 
twice the disk radius of the 'standard' model (180 AU). This model has an inclination of 73° 
and an increased stellar luminosity of 2.9 L0 to fit the observed SED. 
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Finally, we explore the possibility that the disk may be more inclined than ~ 70°. There 
are two ways of increasing the model inclination while preserving the spectral characteristics 
of the source, in particular the high mid-infrared continuum: the first is to decrease the disk 
opening angle and the second is to decrease the mass of the disk. If the disk opening angle 
is decreased by a factor of 2 to 0.09, the inclination can be increased to ~ 76°. However, 
since the disk intercepts less stellar irradiation, the luminosity of the central star has to be 
increased significantly, causing absorbing disk material to increase its temperature to the 
levels of the less inchned disk. Additionally, the disk becomes too flat to fit the near-infrared 
images. If the mass of the disk is decreased to facihtate are higher inclination of 75°, the 
temperature of the disk mid-plane increases to well above 30 K. 



5. Where are the ices toward CRBR 2422.8-3423 located? 

5.1. Ices in the model disk 

The radiative transfer model of CRBR 2422.8-3423 allows the location of the ices to 
be studied in greater detail. In principle, the model cannot directly distinguish between ice 
located in the disk and ice located in foreground material, but significant conclusions can 
be reached by considering the range of realistic ice abundances. Furthermore, differences in 
dust temperatures in the disk and foreground place constraints on the location of the ices. 
Fig. 15 shows the SED for the disk with the foreground material removed. It is seen that 
significant fractions of the ice absorption lines originate in the disk material if the assumption 
of a constant water ice abundance of 9 x 10~^ in the disk and foreground material is correct. 
Quantitatively, the model disk contributes 45% of the observed water ice band for a constant 
ice abundance. 

From Fig. 8 it is seen that the temperature in the model disk nowhere drops below 27 K. 
This effectively prevents any pure CO ice from being present in the model disk, since this 
evaporates at ~ 20K at the densities present in the disk (Sandford & Allamandola 1993). 
A significant change of the disk density structure, such as 1-2 orders of magnitude denser 
material in the mid-plane or a significantly larger outer disk radius, is required in order to 
have pure CO ice be stable in the disk. Even if CO ice is present in a dense mid-plane, it 
is unlikely to contribute to the observed CO ice band. This question can be explored by 
examining the temperature of the dust that can be probed at a wavelength of 4.67 /^.m. Fig. 
16 illustrates the point for the best-fitting model disk. It shows the model optical depth 
of the 4.67 //m CO band as a function of desorption temperature, assuming a constant CO 
ice abundance for temperatures lower than the desorption temperature. The relation has 
been calculated by directly including a component of pure CO inclusions in the dust model 
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for temperatures lower than the desorption temperature. Naturally, if a higher desorption 
temperature is assumed, a deeper CO ice band appears in the spectrum. Therefore, the 
contribution of ice at different temperatures to an ice band centered on 4.67 /xm can be read 
off the relation in the figure. The foreground material is seen as a sharp rise to 20 K followed 
by a plateau, indicating that little absorbing material is present at temperatures between 30 
and 40 K. The contribution from material in the disk takes over only at temperatures above 
40 K. This clearly shows that the model disk can contain no pure CO ice that will show up 
in the spectrum; all ice in the disk observed in the wavelength region around 5 /im must have 
temperatures between 40 and 90 K. This is consistent with the excitation temperature of the 
observed CO rovibrational hnes of 40-60 K as determined by Thi et al. (2002). Inspection 
of Fig. 8 shows that the disk material, which contains ices contributing to the observed 
absorption bands, has densities of 10^ — 10^ cm^^. The disk must be considerably larger in 
order to contribute to the observed pure CO ice band. This is illustrated by plotting the 
CO ice band depth as a function of assumed desorption temperature of the large disk with 
-Rdisk = 180 AU in Fig. 16. It is seen that the large disk shifts the absorbing material toward 
lower temperatures as expected, but only by ~ 5 K, such that the coldest dust in the disk 
contributing to the ice band is 35 K for the large disk. 

If the CO is embedded in other less volatile ice species, such as water ice, the evaporation 
temperature of CO can be raised to 70 K. However, the CO and CO2 ice profiles indicate 
that at most 20% of the CO is mixed with CO2 and another 20% with water. The remaining 
CO ice must reside in the cold foreground dust. The conclusion that none of the pure CO ice 
observed in the spectrum is present in the disk implies that the abundance of pure CO ice in 
the foreground must be very high, 1.1 x 10~^, making CO the most abundant ice species in 
the foreground material. If the CO embedded in water ice is present in the observed average 
abundance of 1.7 x 10~^ (see Table 1) throughout the system, the CO ice abundance in the 
foreground may increase further to 1.27 x 10^^. In order to determine the total depiction of 
gas-phase CO in the foreground, ice species formed by CO, such as CO2 and CH3OH, must 
also be included. Adding the observed average abundance of CO2 of 2.9 x 10~^ and ignoring 
the relatively small contribution from CH3OH, the total amount of CO depleted from the 
gas-phase is as much as 1.6 x 10~^. For a total CO abundance of 2 — 3 x 10~^ w.r.t. H2, this 
leaves at most 1 x 10~^ for gas-phase CO, corresponding to a foreground column density of 
2.2 X lO^^cm"^. This is roughly consistent with the total CO column density (foreground 
+ background) of ~ 5 x lO^^cm"^ derived from millimeter C^^O data (Thi et al. 2002). We 
conclude that the CO gas in the foreground cloud is depleted by a factor of 2-3. 

The assumption that all other ice species evaporate at 90 K is not very realistic and was 
adopted for simplicity. In reality, the ices considered in the model evaporate at a range of 
temperatures between 20 and 100 K, depending on the species and ice composition. Fig. 17 
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shows the variation of the model ice features in the disk if an evaporation temperature of 
45 K is assumed instead. It is seen that the ice bands become ~ 30% weaker for the lower 
evaporation temperature, corresponding to most of the disk ice. The implication is that all 
disk ice features should show significant signatures of heating above 40 K. For CRBR 2422.8- 
3423, this effect may be difficult to observe due to the contribution from cold foreground 
ice. However, as discussed in Sec. 3.3, the red shoulder of the 6.85 /im band indeed provides 
significant evidence for the presence of warm ices along the line of sight. It is important 
to note that most of the disk mass has temperatures below 40 K, but that this material 
is simply not probed by the ice bands due to the high optical depths in the mid-plane 
and the system inclination angle. This illustrates how the location of ices in circumstellar 
disks and probably in most low-mass protostars require a detailed axisymmetric model to 
be adequately constrained. 

In summary, we find that up to 50% of the observed band optical depths of ice species 
that remain in the solid state at temperatures above 40 K, may be due to disk ice. This 
includes water, CO2 ices and CO embedded in water ice. The observed pure CO ice is most 
likely located entirely in the Oph-F core. 
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Fig. 15. — Models showing the effect of inchnation on the ice bands. The models are those 
that best-fit CRBR 2422.8-3423, only with the foreground cloud removed. The inchnations 
shown are from top to bottom: 65°, 69.15° ('standard' model value - see Table 2), 75° 
and 85°. Note the significant changes in fine ratios, especially between the 3.08 //m and 
the 6.0 //m water ice bands. The vertical fines indicate the position of the 3.08 //m water 
stretching mode, the 4.27 fim CO2 stretching mode, the 6.0 //m water bending mode and the 
15.2 /xm CO2 bending mode. 
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Fig. 16.— Optical depth of the 4.67 /xm band of CO for the CRBR 2422.8-3423 model as 
a function of (the assumed) desorption temperature. The abundance is arbitrarily fixed at 
9 X 10~^ whenever the temperature is below the desorption temperature. Both the 'standard' 
model (solid curve) of Table 2 and the model with -Rdisk = 180 AU (dotted curve) are shown. 
The vertical solid and dotted lines indicate the separation between contributions from the 
foreground material and from the disk material for the 90 AU and 180 AU disks, respectively. 
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Fig. 17. — Comparison of a model assuming that the ice evaporation temperature is Tevap — 
90 K (sohd curve) with a model using Tevap = 45K (dashed curve). For CRBR 2422.8-3423, 
a significant amount of the ice has temperatures in excess of 45 K. 
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5.2. The missing water ice librational band 



In Fig. 12, it is evident that the red wing of the 9.7 fim siUcate band is not well-fitted 
by the model. The model red wing is caused by the librational (hindered rotation) band of 
water ice. This band is centered around 12-13 /xm, but the center and width are sensitive to 
ice composition (Hagen et al. 1983). The band strength of the water libration band is almost 
three times higher than the CO2 bending mode, which is located in the same wavelength 
region (Gerakines et al. 1995). Combined with the higher column density of water ice, 
the libration band should be about 10 times stronger than the corresponding 15.2 /xm CO2 
band, although it will only contribute an optical depth of 0.3-0.5 at 12 fj,m due to the large 
width of the band. Laboratory experiments find that the strength of the librational band is 
fairly constant under interstellar conditions (Kitta & Kraetschmer 1983). Careful radiative 
transfer models of other sources compared to, e.g., Spitzer-IRS spectroscopy may show if 
the problem is universal. Recent Spitzer-IRS spectra of the embedded low-mass star HH 46 
IRS, toward which the ice bands are about twice as deep as toward CRBR 2422.8-3423, also 
do not show a water ice libration band as strong as expected (Boogert et al. 2004). The 
hbrational band is sensitive to the dust model, through the grain size distribution as well as 
the shape of the grains. It is possible that a more realistic dust model may provide a better 
fit the librational water ice band. 



6. Observing ices in other highly inclined disks 

Using the model for the CRBR 2422.8-3423 disk, some general predictions can be made 
for observations of ices in other edge-on circumstellar disks. Clearly, the presence of cold 
foreground material complicates any interpretation of ice bands from a disk. However, since 
this problem often cannot be avoided, it is interesting to explore how the ice bands intrinsic 
to the disk may be distinguished from those produced by foreground material. 
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Fig. 18. — Upper panel: Optical depths of the H2O and CO2 ice bands as functions of disk 
inclination for the model of CRBR 2422.8-3423. Lower panel: Line ratios relative to the 
3.08 ^m. water ice band. The horizontal lines indicate the ratios of the extinction coefficients 
of the bands, i.e. the line ratio expected to be observed in a cold foreground cloud. No 
foreground material has been inchided to clarify the effects intrinsic to the disk. The curves 
for the CO2 4.27 /um band have been scaled by a factor of 0.5 (upper panel) and 0.3 (lower 
panel) to fit in the figure. 
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An interesting question to explore is how the ice bands behave with varying inchnation. 
For instance, are disks with an inchnation angle close to 90° more favorable for observing 
ices than disks with lower inclination angles of ~70°? In Fig. 15 the spectra of the best 
fitting model of CRBR 2422.8-3423 viewed at different inclination angles are shown. It is 
apparent that an inclination of 70° in fact produces the deepest ice bands. A disk similar to 
CRBR 2422.8-3423 viewed very close to edge-on shows only weak intrinsic ice features. This 
is an effect produced by the large optical depth through the mid-plane. Disks 1-2 orders of 
magnitude less massive than CRBR 2422.8-3423 will have optical depths approaching unity 
along the mid-plane. Such an edge-on disk will have deep ice bands if ices are present in 
abundance. 

Another observation of considerable interest is that the ratios between different ice 
bands of the same species are not preserved as a function of inclination. To illustrate this 
quantitatively, the model optical depths of H2O and CO2 bands are plotted in Fig. 18 as 
functions of the disk inclination. Common for all ice bands is that the optical depth peaks 
strongly for inclinations between 65 and 75°. For most ice bands a constant plateau is 
reached at higher inclination angles after the main peak. The exception here is the 15.2 /xm 
CO2 bending mode, which drops to a very small absorption depth at the highest inclination. 
Finally, the different bands are seen to peak at somewhat different inclination angles. 

The model optical depths arc also shown relative to the model 3.08 fira water band 
depth in the figure. Clearly, the line ratios of all ice bands relative to the water band change 
dramatically with inclination. Note that this is not an effect of varying ice abundances 
within the disk, since both H2O and CO2 ices are assumed to evaporate at 90 K. The effect 
is especially apparent for the 3.08 //m and 6.0 fj,m water ice bands, the ratio of which drops 
from ~ 10 at an inclination of 65° to ~ 0.7 at an inclination of 74°. At even higher inclination, 
the ratio between the two water ice bands drops to a plateau of 0.2-0.25. 

Note that the varying ratios of the ice absorption bands are, in this purely 2- 

dimensional effect. The 'classical' radiative transfer scenario for ice bands is material in a cold 
slab seen toward a background infrared point source. The important difference for observing 
ices in a disk is that the background source is extended relative to the absorbing material. 
For instance, if the infrared source is due to scattered light, the size of the background source 
may be comparable to the size of the ffarcd disk. There are two mechanisms which produce 
the extended source in the case of a disk. The first is light from the central star and inner 
rim scattered in the surface layer of the disk, and the second is thermal emission from the 
surface layer. Inspection of Fig. 15 illustrates the relative importance of these two effects. 
For the spectrum produced by the highest inclination, a minimum in the flux is observed 
at ~ 8 /im. At shorter wavelengths, scattered light dominates the spectrum, while thermal 
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emission from the disk surface layer dominates at longer wavelengths. This means that the 
15.2 /xm CO2 bending mode is always observed against a thermal source, while the rest of 
the main ice bands arc observed toward a scattered source for a highly inclined disk. At an 
inclination similar to the disk opening angle, the emission from the central star and inner 
rim start to dominate. In this case the entire length of the disk is probed and the deepest 
ice bands are produced. In Fig. 15, it is seen that the inclination where the peak absorption 
of the first three ice bands (the 3.08 //m and 6.0 //m water ice bands and the 4.27 //m CO2 
ice band) occurs, increases with increasing wavelength of the band. This is, however, not 
true for the 15.2 /im CO2 ice band, which peaks at a lower inclination than the 6 /im band. 
This is due to the fact that the first three bands are seen toward a scattering source, while 
the 15.2 /im band is seen toward a thermal source whose geometry is different. 

Obviously, such behaviour significantly complicates the determination of accurate ice 
abundances in disks and other multi-dimensional scenarios. On the other hand, strongly 
deviating ice band ratios seen toward an edge-on disk is a strong signature that the ices 
are intrinsic to the disk. For the case of CRBR 2422.8-3423, no strongly deviant line ratios 
are observed, at least not when compared to other ice sources. In particular, the ratio 
between the 6.0 fim and the 3.08 fim water ice bands is typical of embedded sources. It is 
a long-standing problem that the water ice band 6.0/3.08 ratio is up to two times higher 
than expected. Part of this may be explained by radiative transfer effects analogous to those 
identified here for nearly edge-on disks. However, since most observed icy lines of sight 
exhibit this anomalous ice band ratio, the effect is unlikely to be entirely due to the presence 
of disk geometry. The CO2 bending mode of CRBR 2422.8-3423 can, unfortunately, not be 
compared to the stretching mode at 4.27 /im, which is located outside the spectral coverage 
of Spitzer-IRS. The ratio of the two CO2 ice bands could otherwise constitute an excellent 
tracer of radiative transfer effects in disks. 

The fact that CRBR 2422.8-3423 has an inclination which is exactly optimal for ob- 
serving ices in the disk is an interesting coincidence. Had the inclination been a few degrees 
higher, the mid-infrared continuum flux level would have dropped by an order of magnitude, 
and the ice band ratios would have deviated strongly from the norm. Another important 
conclusion is that pure CO ice is not expected to show up in spectra of disks with radii 
smaller than 100-150 AU and surrounding stars of luminosity > ILq. Any CO ice present 
should show clear signs of being embedded in a less volatile ice matrix, such as a dominance 
of the red water-rich component at 2136 cm~^. 
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7. Conclusions 

We have presented a full mid-infrared spectrum of the inclined disk CRBR 2242.8-3423 
covering all of the strongest absorption features due to interstellar ices, except the A.27 fim 
CO2 band. By fitting an axisymmetric, but effectively 3D, Monte Carlo radiative trans- 
fer model to the spectrum as well as photometric data covering 1.6-3000 /xm and spatially 
resolved near-infrared images, we have been able to constrain the location of ices in the 
disk and in the foreground material along this line-of-sight. The main conclusions can be 
summarized as follows: 

• The use of an axisymmetric continuum radiative transfer model to simultaneously fit a 
high spectral resolution SED as well as high spatial resolution imaging is a powerful tool 
to estimate the physical conditions of dust within circumstellar disks. Combining the 
model with high resolution temperature-dependent opacities further allows a detailed 
study of the icy material present in the disk. 

• The extinction of the near-infrared emission together with submillimeter imaging in- 
dicate that a dense and cold foreground component, containing a significant part of 
the observed ices, is present. Specifically, the line-of-sight toward CRBR 2422.8-3423 
passes through the dense core Oph-F. 

• If a constant ice abundance for all ice species apart from CO is assumed throughout the 

disk and foreground material for temperatures below 90 K, up to 50% of the water and 
CO2 ices reside in the disk. The 6.85 fim feature shows signatures of strong thermal 
processing, indicating that a similar fraction of the carrier, which is possibly NH^, also 
resides in the disk. The shape of the 6.85 /xm band is arguably the strongest evidence 
for the presence of ices in the disk. 

• None of the pure CO ice observed toward CRBR 2422.8-3423 can be located in the 
model disk because the temperature is too high. The presence of a colder self-shadowed 
disk at larger radii is not ruled out, but CO ice in a flat outer disk will not contribute 
to the observed 4.67 /im band. This is because an outer disk with a component along 
the line-of-sight to the central star will also appear in scattered light. According to 
the disk model, only dust with temperatures of 40-90 K and associated gas densities 
of 10^ — 10^ cm~^ can contribute to the observed ice bands. Up to 20% of the CO ice 
is trapped in water ice and can survive in the disk together with an additional 20% 
trapped in CO2 ice. For a constant abundance of these types of CO ice throughout the 
disk and envelope, at most 20% of the total observed amount of solid CO originates in 
the disk. In general, pure CO ice should only be observed in large quantities in disks 
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surrounding low-mass stars if the disks have radii much larger than ~ 150 AU. Indeed, 
Thi et al. (2002) concluded that most of the CO ice is located in the disk, based on 
a larger model disk. Excluding the presence of pure CO ice in the disk required the 
detailed 2D modeling presented here. It should be emphasized that the new Spitzer 
observations indicate that warm rather than cold ices are indeed present in the disk. 

• In general, the high optical depths through the mid-planes of nearly edge-on disks 
prevent cold ices in the disk (< 30 K) at radii of < 200 AU from being directly observed 
in absorption toward a central 1 Lq star. 

• The model does not provide a good fit to the red wing of the 9.7 //m silicate band. 
The discrepancy is due to the absence of the strong librational band of water ice 
centered at 12-13 iim. We are unable to explain the absence of the band since laboratory 
experiments indicate that the strength does not change appreciably under interstellar 
conditions. 

• The infrared source against which ice bands are observed is extended for an edge-on 
disk due to scattering at wavelengths shorter than ~ 8 /im and due to thermal emission 
from the disk surface layer at longer wavelengths. This results in ice band ratios which 
are not equal to the corresponding band ratios of the opacity. The model indicates 
that the ratio between the 3.08 /im and the 6.0 /xm water ice bands as well as the ratio 
between the 4.27;um and the 15.2 fim CO2 ice bands are excellent tracers of radiative 
transfer effects on ice bands intrinsic to an edge-on disk. Ratios between bands of the 
same ice species are potentially powerful probes of the distribution of ices in disks. At 
the same time, anomalous ratios may be used to determine if the ices are located in 
the disk as opposed to a foreground cloud. The opposite is not true, i.e. normal band 
ratios do not indicate that the ices are located outside the disk. 



Support for this work, part of the Spitzer Space Telescope Legacy Science Program, 
was provided by NASA through Contract Numbers 1224608 and 1230780 issued by the 
Jet Propulsion Laboratory, California Institute of Technology under NASA contract 1407. 
KMP was supported by a PhD grant from the Netherlands Research School for Astronomy 
(NOVA). Astrochemistry in Leiden is supported by a SPINOZA graint of the Netherlands 
Organization of Scientific Research (NWO). The authors are grateful to Doug Johnstone for 
providing the SCUBA 850 //m map from the COMPLETE survey and Karl Stapelfeldt for 
comments which improved the manuscript. Finally, we wish to acknowledge the efforts of all 
the people who put many years of hard work into making the Spitzer Space Telescope IRS 
instrument a reality. 



-56- 



REFERENCES 

Aikawa, Y., Umebayashi, T., Nakano, T., & Miyama, S. M. 1997, ApJ, 486, L51+ 
Alexander, R. D., Casali, M. M., Andre, P., Persi, P., & Eiroa, C. 2003, A&A, 401, 613 
Beck, T. L., Prato, L., & Simon, M. 2001, ApJ, 551, 1031 

Bohren, C. F., & Huffman, D. R. 1983, Absorption and scattering of light by small particles 
(New York: Wiley, 1983) 

Boogert, A., Pontoppidan, K., Lahuis, F., J0rgensen, J., Augereau, J.-C, Blake, G. A., 
Brooke, T., Brown, J., DuUemond, C. P., Evans II, N., Geers, V. C., Hogerheijde, 

M., Kessler-Silacci, J., Kncz, C, Morris, P., Noriega- Crespo, A., Schoicr, F., van 
Dishocck. E. F., Allen, L., Harvey, P., Koerncr, D., Mundy, L., Myers, P., Padgett, 

D. , Sargent, A., & Stapelfeldt, K. 2004, ApJ, 154, 359 

Boogert, A. C. A., Hogerheijde, M. R., & Blake, G. A. 2002a, ApJ, 568, 761 

Boogert, A. C. A., Hogerheijde, M. R., Ceccarelli, C, Tielens, A. G. G. M., van Dishoeck, 

E. F., Blake, G. A., Latter, W. B., & Motte, F. 2002b, ApJ, 570, 708 

Brandner, W., Sheppard, S., Zinncckcr, H., Close, L., Iwamuro, F., Krabbc, A., Maihara, 
T., Motohara, K., Padgett, D. L., & Tokunaga, A. 2000, A&A, 364, L13 

Ghiang, E. I., & Goldreich, P. 1997, ApJ, 490, 368 

GoUings, M. P., Dever, J. W., Fraser, H. J., McGoustra, M. R. S., & Wilhams, D. A. 2003, 
ApJ, 583, 1058 

Dartois, E., & d'Hendecourt, L. 2001, A&A, 365, 144 

Dartois, E., d'Hendecourt, L., Thi, W., Pontoppidan, K. M., & van Dishoeck, E. F. 2002, 
A&A, 394, 1057 

Draine, B. T. 2003, ARA&A, 41, 241 

DuUemond, G. P., & Dominik, G. 2004, A&A, 417, 159 

DuUemond, G. P., Dominik, G., & Natta, A. 2001, ApJ, 560, 957 

DuUemond, G. P., & TuroUa, R. 2000, A&A, 360, 1187 

DuUemond, G. P., van Zadclhoff, G. J., & Natta, A. 2002, A&A, 389, 464 



-57- 



Dutrey, A., Guilloteau, S., & Guelin, M. 1997, A&A, 317, L55 

Ehrenfreund, P., Boogert, A. C. A., Gerakines, P. A., Jansen, D. J., Schutte, W. A., Tielens, 
A. G. G. M., & van Dishoeck, E. F. 1996, A&A, 315, L341 

Ehrenfreund, P., Boogert, A. C. A., Gerakines, P. A., Tielens, A. G. G. M., & van Dishoeck, 
E. F. 1997, A&A, 328, 649 

Ehrenfreund, P., Breukers, R., D'Hendecourt, L., & Greenberg, J. M. 1992, A&A, 260, 431 

Ehrenfreund, P., Kerkhof, O., Schutte, W. A., Boogert, A. C. A., Gerakines, P. A., Dartois, 
E., D'Hendecourt, L., Tielens, A. G. G. M., van Dishoeck, E. F., & Whittet, D. C. B. 
1999, A&A, 350, 240 

Evans, N. J., Allen, L. E., Blake, G. A., Boogert, A. C. A., Bourkc, T., Harvey, P. M., 
Kessler, J. E., Koerner, D. W., Lcc, C. W., Mundy, L. G., Myers, P. C., Padgett, 

D. L., Pontoppidan, K., Sargent, A. 1., Stapelfeldt, K. R., van Dishoeck, E. F., Young, 
C. H., & Young, K. E. 2003, PASP, 115, 965 

Falk, M. 1987, J. Chem. Phys, 86, 560 

Gerakines, P. A., Schutte, W. A., Greenberg, J. M., & van Dishoeck, E. F. 1995, A&A, 296, 
810 

Gerakines, P. A., Whittet, D. C. B., Ehrenfreund, P., Boogert, A. C. A., Tielens, A. G. G. M., 
Schutte, W. A., Chiar, J. E., van Dishoeck, E. F., Prusti, T., Helmich, F. P., & de 
Graauw, T. 1999, ApJ, 522, 357 

Hagen, W., Greenberg, J. M., & Tielens, A. G. G. M. 1983, A&AS, 51, 389 

Hogerheijde, M. R., van Langevelde, H. J., Mundy, L. G., Blake, G. A., & van Dishoeck, 

E. F. 1997, ApJ, 490, L99+ 

Jager, C., Mutschke, H., & Henning, T. 1998, A&A, 332, 291 

Keane, J. V., Tielens, A. G. G. M., Boogert, A. C. A., Schutte, W. A., & Whittet, D. C. B. 
2001, A&A, 376, 254 

Kerkhof, O., Schutte, W. A., & Ehrenfreund, P. 1999, A&A, 346, 990 

Kitta, K., & Kraetschmer, W. 1983, A&A, 122, 105 

Leger, A., Gauthier, S., Defourneau, D., & Rouan, D. 1983, A&A, 117, 164 



-58- 



Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 

MuzeroUe, J., Calvet, N., Hartmann, L., & D'Alessio, P. 2003, ApJ, 597, L149 

Natta, A., Prusti, T., Neri, R., Wooden, D., Grinin, V. P., & Mannings, V. 2001, A&A, 371, 
186 

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 

Pontoppidan, K. M., Fraser, H. J., Dartois, E., Thi, W.-F., van Dishoeck, E. F., Boogert, 
A. C. A., d'Hendecourt, L., Tielens, A. G. G. M., & Bisschop, S. E. 2003, A&A, 408, 
981 

Pontoppidan, K. M., van Dishoeck, E. F., & Dartois, E. 2004, A&A, in press 
Sandford, S. A., & AUamandola, L. J. 1993, ApJ, 417, 815 

Schutte, W. A., Boogert, A. C. A., Tielens, A. G. G. M., Whittet, D. C. B., Gerakines, P. A., 
Chiar, J. E., Ehrenfreund, P., Greenberg, J. M., van Dishoeck, E. F., & de Graauw, 
T. 1999, A&A, 343, 966 

Schutte, W. A., Gerakines, P. A., Geballe, T. R., van Dishoeck, E. F., & Greenberg, J. M. 
1996, A&A, 309, 633 

Schutte, W. A., & Khanna, R. K. 2003, A&A, 398, 1049 

Shuping, R. Y., Chiar, J. E., Snow, T. P., & Kerr, T. 2001, ApJ, 547, L161 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 

Smith, R. G., Sellgren, K., & Tokunaga, A. T. 1989, ApJ, 344, 413 

Stapelfcldt, K. R., Menard, F., Watson, A. M., Krist, J. E., Dougados, C., Padgett, D. L., 
& Brandner, W. 2003, ApJ, 589, 410 

Thi, W. F., Pontoppidan, K. M., van Dishoeck, E. F., Dartois, E., & d'Hendecourt, L. 2002, 
A&A, 394, L27 

Thi, W.-F., Zadelhoff, G.-J., & van Dishoeck, E. F. 2004, A&A, in press, astro 

van Zadelhoff, G.-J., van Dishoeck, E. F., Thi, W.-F., & Blake, G. A. 2001, A&A, 377, 566 

Vandenbussche, B., Ehrenfreund, P., Boogert, A. C. A., van Dishoeck, E. F., Schutte, W. A., 
Gerakines, P. A., Chiar, J., Tielens, A. G. G. M., Keane, J., Whittet, D. C. B., 
Breitfellner, M., & Burgdorf, M. 1999, A&A, 346, L57 



-59- 



Warren, S. G. 1984, Appl. Opt., 23, 1026 

Watson, D. M., Kemper, F., Calvet, N., Keller, L. D., Furlan, E., Hartmann, L., Forrest, 
W. J., Chen, C. H., Uchida, K. I., Green, J. D., Sargent, B., Sloan, G. C., Herter, 
T. L., Brandl, B. R., Houck, J. R., Najita, J., D'Alessio, R, Myers, P. C., Barry, 
D. J., Hall, P., & Morris, P. W. 2004, ApJS, 154, 391 

Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 

Werner, M. W., Gallagher, D. B., & Trace, W. R. 2004, Advances in Space Research, 34, 600 
Wood, K., Wolff, M. J., Bjorkman, J. E., & Whitney, B. 2002, ApJ, 564, 887 



This preprint was prepared with the A AS IM^jX macros v5.2. 



